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Abstract 

In this paper we analyze a two-dimensional discrete model of nearest- 
neighbour Lennard-Jones interactions under the microscopical constraint 
that points on a lattice triangle maintain their order. This can be un¬ 
derstood as a microscopical non-interpenetration constraint and amounts 
to the positiveness of the determinant of the gradient of the piecewise- 
afhne interpolations of the discrete displacement. Under such a constraint 
we examine the continuum fracture energy deriving from a discrete-to- 
continuum analysis at a scaling where surface energy is preponderant. 
We give a lower bound by an anisotropic Griffith energy. This bound 
is optimal if the macroscopic displacement satisfies some opening-crack 
conditions on the fracture site. We show that if such conditions are not 
satisfied then the computation of the energy due to continuum cracks may 
involve non-local effects necessary to bypass the positive-determinant con¬ 
straint on crack surfaces and at points where more cracks meet. Even 
when the limit fracture energy may be described by a surface energy den¬ 
sity, this may depend on the crack orientation both in the reference and 
in the deformed configuration. While these effects lead to very interesting 
analytical issues, they call into question the necessity of the determinant 
constraint for fracture problems. 

Keywords: Lennard-Jones potentials, variational theory of Fracture, dis- 
crete-to-continuum analysis, T-convergence. 


Introduction 

In this paper we analyze on one hand continuum descriptions for fracture 
energies arising from discrete systems of particles linked by interatomic 
interactions, and on the other hand computational problems deriving from 
microscopic constraints. In the simplest model of atomistic interactions, 
the behavior of a collection of N particles is governed by an energy that 
can be written as a sum of pair contributions; namely, it can be set in the 
form 

En{{u.})= Y. ( 0 . 1 ) 

where Ui is the position of the i-th atom, ||Mi — Uj\\ the distance between 
the corresponding pair of atoms, and J is an interatomic potential, which 
is strongly repulsive at short distances and mildly attractive at long dis¬ 
tances. The most common choice for J is the Lennard-Jones potential 

JA) = ^ - 5 ( 0 . 2 ) 
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(ci,C2 > 0). Note that in principle the total energy of the system is 
bounded from below by N'^ min J; however, a more refined estimate shows 
that it scales exactly as the number N of atoms. This is in agreement with 
the intuition that gronnd states for En arrange approximately in a regular 
lattice as N increases {crystallization), so that the energy contribution 
of each particle is essentially described by the interaction only with its 
neighbours. At the same time it suggests that models of crystalline solids 
can be derived directly from such atomic interactions. 

The crystallization for Lennard-Jones interactions in the general con¬ 
text describe above is still an open problem. An important contribution 
has been recently given by Theil ]22\, who has studied a slightly weaker 
version of crystallization, proving in two dimensions the optimality under 
compact perturbations of the “minimal” triangular lattice; i.e., the one 
for whose lattice spacing f it is minimal the average energy per particle 

eir)= ■^(^11*11)’ (0-3) 

fceT\{o} 

which describes the energy of a single particle in the lattice rT due to the 
interaction with the other particles, T being the unit triangular lattice. 
Once crystallization is achieved, another important issue is whether it 
is maintained on states other than ground states; i.e., whether to (small) 
macroscopic deformations there corresponds uniform displacements at the 
atomic level {Cauchy-Born rule-, see Friesecke and Theil |21p. 

In order to examine the behavior of atomistic systems far from ground 
states, under the hypothesis of crystallization we may consider the energy 
related to a density p. Noting that p is proportional to r~", where n is 
the space dimension (usually, n = 2, 3), this can be expressed as 

f{p) = p(Cip^^/" - 


where Ci, C 2 are determined by computing the energy of a single particle 
in a uniform lattice (in the case of a triangular lattice f{p) = pe{r), where 


e is defined in (0.31 and r is the lattice spacing corresponding to the density 


p). This function / is non-convex in an interval (0,po), which highlights 
a phase transition at low densities, and suggests that large deformations 
involve a change in the crystalline structure which is instead achieved close 
to ground states. 

From the standpoint of Continuum Mechanics, a description with an 
hyperelastic bulk energy is expected to hold close to ground states, while 
the same is expected to fail for a class of large deformations. In the one¬ 
dimensional case this can be achieved by introducing a fracture energy. 
Following Truskinovsky [24], in this case, given N particle positions Ui, 
these can be parameterized with i = 1,... in such a way that Ui > Ui-i, 
and write the energy as 

i?iv({wi}) = ^ Ji'^3 — “0- 

j>i 

As N increases ground states tend to arrange regularly on a lattice that 
we may suppose to be Z; i.e., we may suppose that the one-dimensional 
energy per particle 

e{r)^ 

fcGN, fc>0 

has its minimum at r = 1. In order to introduce a macroscopic deforma¬ 
tion gradient, we can now scale and re-parameterize the same particles 
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on eZ, where e = ^; i.e., Ui = -u{ie)^ so that they all can be seen as 
discretizations of functions defined on a single interval [0,1]. If we scale 
the energy as 

j>i 

then the ground states are discretizations of the identity on [0,1]. If we 
let N increase (or e —^ 0), we may highlight two regimes: 

{bulk scaling) under the hypothesis of small perturbations u{x) = x + 
Sv{x) with small 5,we have 


u{{i + k)e) — u{ie) 
e 


k + k5v'{ie), 


so that, Taylor expanding e at its minimizer 1, 


i5iv({ui}) = Aimine + 


|w'|^dx + o((5^). 


{surface scaling) if the macroscopic u is discontinuous at a point pa¬ 
rameterized by i then we have an increase from the ground state energy 
of the order J(-|-oo) — mine, which gives a Griffith fracture energy at each 
such point. 

This argument can be made rigorous, and it gives: 

• eiastic behavior ciose to ground states, with a iinearized description 
given by the linearization of J at ground states; 

• a brittle fracture energy depending on the depth of the well of J 
with respect to the infinity; 

• opening fracture: the possibility of a parameterization with increas¬ 
ing u implies that fracture may only open up, providing a natural non¬ 
interpenetration condition; 

• surface relaxation: on external boundaries and on internal fracture 
sites the asymmetry of the atomic arrangements gives an additional sur¬ 
face term, highlighting a microscopic rearrangement close to those sur¬ 
faces. 

We note that all the features above can be included in a single analysis 
by choosing S = in the notation above, so that (up to additive 

constants) bulk and surface terms have the same scaling. This has been 
done by Braides, Lew and Ortiz [9]. We also remark that an analysis of 
local minima and of opening cracks suggests a cohesive fracture energy 
density of Barenblatt type, which is not directly given by the analyses 
above. However some cohesive theories can be shown to be “equivalent” 
to the ones obtained by the limit analysis (see Braides and Truskinovsky 

m and 0). 

The same study for the two or three-dimensional case presents greater 
difficulties, mainly because a natural parameterization with increasing 
functions is no longer possible. Hence, some simplifications have been 
made for this model with the scope of maintaining the features of the com¬ 
plete system of interactions and at the same time allow for an analytical 
study. A first assumption is to consider displacements as a perturbation 
from a ground state for which crystallization holds; i.e., in dimension two 
a perturbation from a state parameterized on the triangular lattice T, or 
rather on a bounded portion A of T (this corresponds to take N as the 
numbers of points in A in the previous notation). For such perturbation 
it makes sense to assume that only short-range interactions be taken into 
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account. The range of such interaction can be indexed with a subset S of 
T \ {0}. The energies replacing En can be written as 

fa{u) = y ] y ] Jk{\\ui — ujII), 

kGS i,jGA,i—j=k 

where Jk represents the energy between points whose parameters i,j dif¬ 
fer by k in the reference lattice. The simplest choice is considering only 
nearest-neighbour interactions on T, with S as the unit vectors of the 
triangular lattice, and Jk independent of fc G S' a Lennard-Jones potential 
with minimum in r. In this simplified model interactions are minimized 
when the corresponding Ui — Uj are at distance r, thus recovering uniform 
deformations. Unfortunately, such a simplihed model admits many ad¬ 
ditional minimizers, as all deformations which are piecewise homotheties 
with ratio ±r. In fact, if we compose a uniform deformation (e.g. a homo- 
thety) of ratio r with “folding” along a line of points in T, the resulting 
nearest neighbours still are at the “minimal” distance r. 

In order to prevent undesired “foldings” at a discrete level without 
considering longer-range interactions, Friesecke and Theil [3T] proposed 
to add a three-point condition on neighbouring nodes. In the case of a 
triangular lattice, this condition simply amounts to requiring that 

det Vu > 0, 

where u is the affine interpolation of the function defined on the vertices of 
each triangle. In this way only discretizations of homotheties with positive 
ratio equal to r are minimizers of the energy. In their paper, Friesecke and 
Theil treat the seemingly “un-natural” case of energies parameterized on a 
square lattice. Actually, when finite-range lattice energies are considered, 
the choice of the underlying lattice is a matter of convenience, and the 
square lattice is the simplest where to consider at the same time nearest 
neighbours and next-to-nearest neighbours to highlight the possibility of 
non-uniform ground states {failure of the Cauchy-Born rule). 

In this paper we treat a two-dimensional system of nearest-neighbour 
Lennard-Jones interactions with the positive-determinant constraint fo¬ 
cusing on the surface scaling. We will use the terminology and techniques 
of F-convergence applied to a discrete-to-continuum analysis [5]. In this 
framework we examine the overall behavior of the energies Fa as the size 
of A diverges, by considering A = n T), with a fixed bounded 
domain in and using U n eT as the set of parameters. The scaled 
energies we are going to examine will be of the form 

F,W= '•'( 11 ^ 1 ). 

i,j£Q,n£T,\i—j\=E 

where Ui is the value of the discrete function u at the node ei € U D eT, 
and the piecewise-afEne interpolation of u on the triangulation related to 
eT is supposed to satisfy the positive-determinant constraint. The scaling 
e heuristically can be explained by considering, as in the one-dimensional 
case, the contribution of a set of indices I where ||ui —Uj|| » 1 and noting 
that under this scaling the finiteness of the energy asymptotically implies 
that they have the dimension of a line, so that they can be regarded as 
an interface. Under these assumption we will address the two issues 

• determine whether some condition of “opening crack” still hold in 
the two-dimensional case; 
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• characterize a limit continuum surface energy defined on functions 
defined on ff. 

The other two issues present in the one-dimensional analysis; i.e., the 
characterization of the bulk energy close to ground states and surface 
relaxation have been separately addressed in slightly different hypotheses 
by Braides, Solci and Vitali [12] (for the bulk analysis) and Theil [23] (for 
the external surface relaxation). An analysis of the small-displacement 
regime has been carried over by Friedrich and Schmidt [1711181 [19] . 



Figure 1: A reference configuration with fracture site and its macroscopic nor¬ 
mals, and its underlying triangulation 


In order to makes the analysis clearer we will scale the energies by an 
additive constant as 

Fs{u) = - ^1^—minjj, 

i,jGnneT,|i—j|=£ 

so that the energy density is always positive. As a first remark we will note 
that, using a result by Chambolle, Giacomini and Ponsiglione [14], gra¬ 
dients of limits of sequences (ue) with equi-bounded energy are piecewise 
rotations with an underlying partition of Q into sets of finite perimeter. 
This allows us to focus on the boundaries of such sets, where we have 
a normal u on whose two sides we have the values u^{x) of u and two 
rotations among those labeling the sets of the partition. Note that 
from the standpoint of the microscopical triangulations the interfaces are 
the limits of triangles of side-length e which are deformed by into tri¬ 
angles with one side (actually two) of diverging length but with the same 
ordering of the vertices (this corresponds to the positive-determinant con¬ 
straint). If only a layer of triangles is deformed that gives a limit interface, 
then this gives a relation between u, u^{x) and In the simplest case 
when u is orthogonal to one of the unit vectors of T (we call those vectors 
coordinate normals), this relation reads as 

(u'''(x) — u~(x}, v) > 0 . 

This can be regarded as an opening-fracture condition in the “finite” 
case: note that in the “infinitesimal case” when u is a small variation of 
a ground state; i.e., u = id + 5v with 5 « 1, then this condition reduces 
to the usual infinitesimal crack-opening condition 

(u''"(®) — v~ {x), v) >Q 
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Figure 2: A deformed microscopical configuration satisfying the determinant 
constraint, and its macroscopic limit with the corresponding rotations 


since = id + o(l). If v is not a coordinate normal then the opening- 
crack condition gets more complicated, due to microscopic anisotropies 
of the lattice, which disfavor cracks not orthogonal to lattice directions. 
The same anisotropies appear, as usual for lattice systems, in the form 
of the surface energy density. However, the situation described above is 
not the only possibility, since more than one layer of triangles may be 
“strongly deformed”. This gives a higher energy on the interface, but 
relaxes the constraints on j/, u^{x) and . Moreover, additional energy 
contributions may be given by points where three or more interfaces meet; 
in this case, even though the opening-fracture condition above may be 
satisfied on each interface, the system of interfaces may be incompatible 
with the positive-determinant condition in their common point. It must 
be observed that by introducing an high number of extra interfaces at 
the discrete level, the finite opening-fracture condition can be removed 
altogether, at the expenses of a complex non-local form of the Hnal energy. 

We Hnally note that, even when non-local effects are neglected, the 
fracture energy density depends not only on the crack orientation v and 
the crack opening — u~, but also on the image of the crack in the 
deformed configuration, described by the tangential derivatives of on 
both sides of the crack. Such types of energies seem to be of interest 
in themselves, and the corresponding analytical techniques still relatively 
undeveloped (see my 


1 Some notation and preliminary results 

In the sequel we denote by (■, ■) the scalar product in and with | • | the 
usual euclidean norm, both for vectors in R^ and matrices in A1(2 x 2; R), 
the meaning being clarified by the context. Moreover we will also use the 
symbol {A, y) to denote the duality pair for A £ Af(2 x 2; R) and j/ £ R^. 
For every x,y £ R^, [x,y\ denotes the segment joining x and y. 

The functional space involved in our discrete-to-continuous analysis 
will consist of measurable vector-valued functions whose components are 
special functions with bounded variation. Such functions are commonly 
used in the variational theories of Fracture (see e.g. the seminal paper by 
Francfort and Marigo [IS], the book by Bourdin, Francfort and Marigo [5] 
and the review paper my 
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Definition 1.1. Given fl C a fixed bounded open set, a function 
u £ is a function with bounded variation in Q, denoted by 

u £ BV (fi, R^), if its distributional derivative Du is a vector-valued Radon 
measure on Q,. 

We say that u is a special function with bounded variation in fl, 
and we write u £ SBV(fl,R^), if u £ and its distributional 

derivative Du can be represented on any Borel set A (Ztd as 

Du{A) = / Vu{x)dx-\- / {u'^{x) — u~(x)) iSi Uu{x) d'H^{x). 

JA JAnSu 

for a countably 'Hf -rectifiable set Su in that coincides 'H} -almost every¬ 
where with the complement in Q of the Lebesgue points of u. Moreover, 
\7u{x) is the approximate gradient of u at x, Vuix) is a unit normal to 
Su, defined for 'Hf -almost every x, and u~^(x),u~ (x) are the traces of u 
on both sides of Su (according to the choice of Uu(x)). Here the sym¬ 
bol ® stands for the tensorial product of vectors; i.e., for any a,b £ 

(a (g) ;= Uibj. 

For fine properties of BV and SBV functions and a rigorous definition 
of all the quantities introduced above we refer to [1], Ch. 4. SBV functions 
enjoy a good compactness property that here is stated in its simplest form, 
suited for our functionals. 

Theorem 1.2. Let (u„) C SBV(fl,R^) be satisfying 

supjy |Vu„|’’dx + ||w„||j;,oo(n,K 2 ) + <+ 00 , (1.1) 

for some p > 1. Then there exist a subsequence Uu(k) o^nd a function u £ 
SBV{Q,R'^) such that Uu (fe) -A u in measure. In addition D°‘Uu;k) D°‘u 
and D^u„[k) D^u weakly* as measures. 

As a consequence we easily get that the limit function u in the previ¬ 
ous theorem has the further property that /_ |Vu„|^ dx-|- HuHnoom + 
n\Su) < + 00 . 

Moreover estimate allows to discuss separately the lower semi¬ 

continuity of suitable bulk and surface integral involving, respectively, 
the absolutely continuous and jump part of the generalized distributional 
derivative. 

In order to treat lower semicontinuous bulk integrals in the vectorial 
setting it is necessary to recall the notion of quasiconvexity and quasicon- 
vex envelopes. 

Definition 1.3. We say that a lower-semicontinuous integrand ip : R^^'^ 

R U {±C)o} is quasiconvex if 

ip{A)C^{yi)< f Ip {A-£ Vip{x)) dx 

Jn 

for every A £ R^^^ and ip £ Wq’°° (fl,R^). In addition, for any f : 
R^x2 u {±oo}, the quasiconvex envelope of f is defined as 

P'^{A) := sup{(ji(A) : g quasiconvex, <?</}. 

In case / is real valued and continuous one can give a precise descrip¬ 
tion of its quasiconvex envelope. 
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Proposition 1.4. Let / : —>■ R be continuous. Then 



forAe M(2 X 2;R). 

The following theorem deals with lower semicontinuity along sequences 
of SBV functions weakly converging [4]. 

Theorem 1.5. Let (m„) C SBV(fl,R'^) be satisfying 

supjy^ |Vun|^ dx +< +00 (1.3) 

for p >1. If Un ^ u in measure, then u € SBV{Q,,R^). Moreover, for 
every quasiconvex integrand ip : R^^^ [0, +oo) such that 

\iPix)\<c{i + \xn 

for every X € R^^^ and C a positive constant, and for any norm p inR?, 
there holds 

f Ip (Xu) dx < liminf [ ip{Xun) dx, (1.4) 

Jn " Jn 



< lim inf 

n 





The final tool we present is a ‘slicing result’ that allows to treat 2- 
dimensional energies by reducing to 1-dimensional ones [1]. 

Before entering into details we introduce some notation. For ^ 
let ~ {y € R^ : (y, = 0} be the line through the origin orthogonal 

to II y G and i? C R^ we set 

-{tGR-.y + t^GE}. (1.5) 

Moreover, if u : i? —>■ R^ we define the function —>• R^ by 

M^’“(t) :=-u(y + tC)- (1-6) 

Theorem 1.6. (a) Let u G SBV{Ll,R^). Then, for all G the func¬ 
tion belongs to SBy(fl^’^,R^) for -almost every y G D^. Moreover 
for such y we have 

u^'^{t) = {Xu{y + tC)i C) for almost every t G , 

S{u^’^) = {tGR:y + t^G S„}, 

M^’^(t±) = m’^(i/+ t^) or {t±) = {y-\-1£,), 

according to the cases {v-ayQ > 0 or {i’u,£,) < 0 (the case {yu,Q = 0 being 
negligible) and for all Borel functions g 

g{t)d'H^{y) = f gix)\{yu,0\d'H'-. 

tes(,u(,v) 

(b) Conversely, if u G L°°{Q,R^) and, for all ( G { 61 , 62 } and for 
TL^-almost every y € 11^, G SBV{Ll^''‘'jR^) with 

[ ([ + #(S'(w^’^))) d'H^(i/) <+ 00 , 

then u ^ 



We introduce now a generalization of a rigidity type result in the SBV 
setting due to Chambolle, Giacomini, Ponsiglione (see [11]). Before going 
through the statement we need some preliminary definitions. First we 
recall the notion of Caccioppoli partition of a domain (see e.g. ESI). 
Definition 1.7. We say that a collection of pairwise disjoint measurable 
sets {Eh}hen is a Caccioppoli partition of fl ifVJhenEh = ff and 

PeT{Eh) < + 00 . 

h£TI 

Moreover, given any rectifiable set K d Pi we say that a Caccioppoli 
partition of Q is subordinated to K \ Uhend*Eh) = 0. 

The following result is the 2-dimensional version of the one contained 
in [T4| . 

Theorem 1.8. Let u € with 'H^{Su) < +oo and \/u{x) € 

50(2) for almost every x & PI. Then there exists a Caccioppoli partition 
{Eh}hen of PI subordinated to Su such that for almost every x G PI 

u{x) = RhX + qh on Eh 

where Rh G 50(2) and qh G . 

2 Formulation of the problem 

In the sequel PI will be any open bounded open set in R^. With Hxed 
discretization step e > 0 the reference lattice is given by := £{ri^, 
where rfi = (1,0), = (1/2, \/3/2). We introduce also the following 

notations 

rf — rj'^ — rj^, S = ± 77 ^, ± 77 ^}. 

Note that S is the set of unitary vectors in the lattice Li and for each 
± G i + S is the set of its nearest neighbours. 

We define also the set D of coordinate directions as 

D = {77^ : 77 GS}. 

Dropping the dependence on e whenever no confusion may arise we 
will use the symbol T to denote any triangle with vertices in and sides 
of length E and 7/ will denote the sets of all such triangles. As already 
pointed out in the Introduction the choice of the lattice relies on the fact 
that Li is the simplest Bravais lattice in dimension 2 compatible with a 
Cauchy-Born hypothesis. Before introducing the precise dehnition of the 
functionals object of our analysis we list here some notation used in the 
following sections for sets of indices or triangles: 

Me{Pl) = {(i, j) G L^x : [i, j] C n, |i - j| = £, i + j}, 

Tfi{P) = {T G% : T CPl}, 

where the symbol i + j stands for the standard lexicographic order in R^. 
Thanks to this choice any pair of nearest-neighbouring indices is counted 
only once in the energy contribution. 

With this notation given a discrete vector-valued displacement u : 
!/£ n —7 R^ we consider its associated energy 

(i,j)eA7e(n) ^ 
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where J : [0, +oo) —> [0, +oo) is a continuous function satisfying the 
following structure properties: 

(i) min J = J{zo) = 0; 

(ii) lim J{z) = Joo > 0; 

z—f + oo 

(iii) lim J(z) = +oo; 

(iv) for any 5 > 03 c« such that Jiz) > Ci(z — zq)’^ for \z\ < 5. 

In what follows for the sake of computational simplicity we simply 
assume zq = 1. Clearly, up to a scaling argument, the analysis remains 
valid for the general case. 

As customary, in order to pass from discrete systems to a continuum 
formulation, it is convenient to identify a function {u(i)}igLjnn C 
with an element of 

Definition 2.1 (interpolation). Given a discrete vector-valued function 
M : n if —>■ we define its interpolation in the whole Q, as a function 

coinciding on each triangle T G with the linear interpolation of the 

values in its vertices. 

Note that an interpolation is not uniquely defined on triangles close 
to do.. This will not affect our arguments since convergences are always 
considered in the interior of SI. We will identify a discrete function with 
its interpolation and maintain the same notation for both the discrete and 
continuous version, the notation being clarified by the context. 

Remark 2.2. In this setting another common procedure is to identify 
a discrete function with a piecewise-constant element in I/^(SI,R^) by as¬ 
signing constant value u(i) on the cell {i 3- -|- pLrf) : A, p G [0,1]}, 

i G n SI. In fact, none of the results stated in the following would 
be affected by this choice (see the discussion contained in m)- Consid¬ 
ering discrete functions as continuous piecewise-afhne ones allows us to 
formulate the orientation preserving constraint in terms of the standard 
determinant of Vu. 


We are now able to introduce the class of admissible functions; i.e., a 
proper class of vector-valued functions on the lattice Lg that are orienta¬ 
tion preserving: 

Ade{id) = {u : Le n n —^ R^ : det(VM) > 0 a.e. in Tfi{Q)}. 

Finally, we define the functional on L^(SI,R^) as 


Ez{u) 


< (i.j)eArE(n) 

[ 3-00 


if u G Ads(fl) 

if u G L^(n,R^)\Ad,(Q). 

( 2 . 1 ) 


Remark 2.3. Note that energies defined in (2.11 account only for inter¬ 
actions well contained in SI and the orientation-preserving constraint is 
not imposed a priori on an affine extension of u on the triangles inter¬ 
secting dfl. Actually, if one is interested in minimum problems endowed 
with boundary conditions (together with some perturbation or fidelity 
terms), a standard procedure in the framework of F-convergence is to fo¬ 
cus on the ‘principal’ part of the total energy, neglecting at first boundary 
data. Once this task is accomplished one can further deal with the initial 
problems up to modifying the limit energy according to the contribution 
arising from recovery sequences with correct boundary datum. 
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To proceed further with our analysis we need to fix a convergence on 

L\n,R'^). 

Definition 2.4 (discrete-to-continuum convergence). According to the 
our identification of discrete functions with interpolations, given Ue,u £ 
we say that —>■ u, and we write simply Wg —>■ u, if we have 
suPg ||wg||ioo(Q_R 2 ) < +00 and linig^o+ ll“£ “ ^^ILi{n,R2) = 0. 


3 Compactness and a lower bound 


In this section we first provide a description of the domain of any F-limit 
functional of the discrete energies defined in (2.11. As already mentioned 
in the Introduction this space consists of piecewise rigid deformations u 
with finite crack energy in the sense of Griffith’s theory; i.e., < 

+ 00 . 


As a second result we exploit geometric measure theory techniques to 
establish a lower-bound estimate for the T-liin infg (u) without imposing 
any a-priori hypothesis on the deformation u. Note that in the next section 
this bound will be proved to be optimal for the class of piecewise rigid 
deformations ‘with opening fracture’ in the sense of the anisotropies of the 
reference lattice. For such deformations the limit fracture energy is simply 
governed by an anisotropic GrifRth-type energy density which reflects the 
anisotropies of the underlying triangular lattice. 


Proposition 3.1. Let {ug} C Ade{Q) be such that liminf i5g(Mg) < +oo 
and 

SUp||Mg||_L=o(fl,K2) < +00. (3.1) 

e 

Then there exists a Borel function u such that, up to subsequences, Us ^ u 
in Z/^(fl,R^). Moreover, 

(i) (finite Griffith fracture energy) u £ 5i3F(fl,R^) with TC^^Su) < 
+oo; 

(ii) (piecewise rigidity) there exists a Caccioppoli partition {Eh}heN sub¬ 
ordinated to S„ such that for almost every x £ Eh we have u{x) = 
RhX + qh for suitable Rh € 50(2) and qh € R^. 


Proof. As a first step we observe that any pair (i, j) £ A4(f^) belongs 
to two different triangles having the segment [i, j] as a side. Hence if we 
take into account a factor 1/2 we may estimate the energies E^(u) from 
below with the integral functionals obtained as a superposition of gradient 
energies indexed by the triangles T varying in Actually, for any 

B open set compactly supported in G and for any v £ Adg)!!), we have 


E,{v) > 


> 


e 

W\ 


7 _ 1 '' T 


Ter/(n) fc=i 
3 


•\/3e 


^,7(KVu,77")|)dx. 


(3.2) 


According to this standpoint, we are led to considering integral functionals 
with energy density J{A) ;= Hence, with fixed any 

s £ (0,1) we distinguish ‘good’ or ‘bad’ triangles T depending whether 
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the triangle energy overcomes the threshold sJoo or not. Using a standard 
separation of scales (see [8]) we set 


T, 


>sJ„ 


and define £ SBV(fl,R^) as 




Ue if T€Z\ X, 

I if T £1, 


(3.3) 


with 7 denoting the identity deformation. By construction the sequence 
(ug) lies in SBV(fl,R^) and for any e > 0 its jump set Sv^ is contained 
in the boundary of the union of the triangles in Ze- 

Moreover, for any open set B compactly supported in inequality 
(3.21 can be rewritten in terms of as 


JiVVe) dx +CsJaoe#{Ze) (3.4) 

v3£ Jb 


for a positive constant c. This implies at once that the functions and 
Ve differ in a set of vanishing Lebesgue measure so that it is enough to 
prove that is compact in SBV{Q,,R'^). 

Thanks to hypothesis (iii) on J the set 

IC = {A£ M(2 X 2;R) ; J(A) < s J^} 


is a compact set in M(2 x 2;R) and this easily provides the estimate 


sup||Vne||i2(B,R2) < +00. 


On the other hand hypothesis (iv) ensures that there exists a constant 
c = c(s) such that (3.41 can be further sharpened as 

1 f ^ 

E,{ue)>- Vc(s)(|(Vu„? 7 '')| - 1)Mx + csJoc 77^(5„, nS). (3.5) 

^ “'s 

Arguing by polar decomposition we also deduce that for any A £ 1C there 

exists R = R{A) £ 50(2) such that 5I)fc=i(K^i'7*^)1 “ 1)^ ^ 

thus 


E. 


(w^) > — / dist^(Vu^,50(2))dx + csJoo'H^(5„, (hB). (3.6) 

£ J B 


Taki ng in to account hypothesis (3.11 a straightforward application of The- 
and the convergence to 0 of Ug — Me, yields that there exists 


1.2 


M £ SBV{Q,R'^) with 'R}{Su) < +00 such that Mg —^ m in L^(S 7 ,R^). 

To prove (ii) we take advantage of a relaxation argument together 
with some rigidity estimates. Indeed, as a first step we prove that \/u{x) £ 
SO{2) for almost every x £ Q. To this end denote V'(A) := dist^(A, 50(2)) 
and estimate the right-hand side of (|3.6[) as 


eE,{u,) > c(s) f i;‘^‘^{Vv,)dx + cesJ^U\S^^nB). (3.7) 
J B 

Passing to the liminf and using (1.41 we get at once that 


L 


i/’^'^(Vm) dx = 0. 
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Hence {Vu{x)) = 0 almost everywhere in 11 and by applying Lemma 


3.3 with A = \/u{x) we also get that the approximate gradient Vw is a 


rotation for almost every x £ B. Once this fact is established (ii) follows 
by applying Theorem |1.8| recursively to u with 12 replaced by any open 
set B compactly supported in Q. and then letting B invading 12. □ 


Remark 3.2. Hypothesis ( 3.1 1 is essential to deduce a compactness result 
for sequences equibounded in energy and avoids the un-physical situation 
of particles escaping to infinity (modelled by arbitrarily large transla¬ 
tions). 


The following lemma completes the proof of Proposition |3.1| This 
kind of results are widely used in variational problems involving crystal 
microstructures. 

Lemma 3.3. Let 4’{A) := dist^(j4, S'0(2)) and let A £ M{n x n, R) be 
such that 'ip'^'^{A) = 0. Then A G SO(2). 


Proof. For A G A1(2 x 2,R) fixed, taking into account Proposition 1.4 
let {(pk)k contained in TTq’°°( f2, R^) be such that 

lim j- dist^(j4 + \/ifik{x), SO{2)) dx = ip‘‘'^{A). 

^ Jn 

Setting for almost every x G Bk{x) = argmindist(A-|-V</9fc(a;), 50(2)), 
we have 


[ \V^kfdx<2[ \A-Bk{x)\^ dx+2 f dist 2 ( 4 -bV(^fc( 2 ;), 50 ( 2 ))d: 
Jn Jn Jn 


X < c. 


Hence it is not restrictive to assume that i/jj, ^ i/Poo and B^ B^o weakly 
in ITq’^( r2, R^) and L^(f2, R^), respectively. By lower semicontinuity, using 
that %j>’^'^{A) = 0, we also get 

0 = •i/;'^‘'(T) = lim / \A+Vifik{x)—Bk{x)\^dx>-[- |4-|-Vv5oc.(a:)—Boo(a;)|^ dx. 
Jn Jn 

Thus V' (j4 -I- V</3tx>(a:)) = 0 for almost every x G LI and the claim follows 
from the classical rigidity theorem in TT^’^(f2,R^) and the fact that ip G 
Wo^’^(f2,R2). □ 


We underline that in the proof of Proposition |3.1| any choice of the 
parameter s G (0,1) ensures the validity of estimate ( |3.7[ ). Moreover, 
assuming that converges to a given u in L^(f2;R^), any sequence = 
Ue(s), defined in (3.3|, still converges to u weakly in SBV{Ll,R.^) and the 
energy contribution is proportional to Since Sv^. is a sequence of 

rectifiable sets with normal coordinates this suggests that in the passage 
to the liminf on e on the surface part of the right-hand side of (3.71 
we may obtain a lower bound with a surface-type energy maintaining 
the simmetries of the hexagonal lattice. In the sequel we built the right 
surface energy arguing by pairwise interactions along lattice directions 
and exploiting more refined techniques in geometric measure theory. 


Before entering into details we need to introduce some more tools. Let 
i/) : R^ —>■ [0, -boo) be the 1-homogeneous map such that {ip < 1} coincides 
with the convex hull of the set S of the unitary vectors of the lattice Li 
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and let t/i* : —>■ [0, +oo) its dual norm; i.e., ip* is the polar function of 


i/) defined as 


sup 

151=1 wiO 


An easy computation shows that 

ip*{iy)= sup 
fc=1.2,3 


and tp* is the 1-homogeneous functions whose unitary ball is the convex 
hull of the coordinate directions D scaled by a factor 2/\/3- In addition 
the following lemma holds true. 

3 

Lemma 3.4. If ip* is as above then 2ip*{i') — \ {h',rjk)\ for all G R^. 

k = l 


Proof. A direct computation shows that the inequality holds true as an 
equality for u G D. Hence one can argue locally in each sector of amplitude 
7r/3 by using the linearity of ip* on such portions of R^ and the result for 
the coordinate directions. The claim follows from the invariance of the 
inequality under rotations with angle tt/S. □ 


By means of Lemma [3.4| in the next proposition we will provide a lower 
bound on T-li minf of Eg by a suitable anisotropic surface energy. Since 
by Proposition 3.1 the P-liminf of Eg is finite only on SBV(Q,M.^) we 


prove the estimate in that functional space. 


Proposition 3.5. 

it holds 


Let (p = Joo(4/%/3)V’*, then for any u G SBV{Ll,S.‘^) 

P-liminf ife(M) > [ ip{uu) dH^. (3.8) 

2^0+ Js^ 


Proof. Let {ug}g C L^(fl;R^) be fixed with Ug —>■ u in L^(fl,R^) and 
l|w£||i°°(n,K 2 ) + Eg(ug) < c. It is not restrictive to assume also that 
liminfg_>g+ Eg{ug) = limj_^Q+ Eg{ug) and Ug ^ u almost everywhere in 
If. Note that by Proposition |3.1| u is a piecewise rigid deformation with 
H^iSu) < +00. 

We will proceed by a slicing technique, performed only in the lattice 
directions rji,ri 2 ,rj 3 , in order to obtain the estimate 


lira Eg{ug)>^^ f \{u„{x),r]k)\d'H\x). (3.9) 

^-0+ ^ 


To this end we observe that we can split the energies Eg by accounting 
separately the contribution of pairs (i, j) with j — i = ErjkS and letting 
k varying in {1, 2, 3}. Hence, setting for v G Adg(fl) 

Eg\v)= Y. + 

iGl/enn 


we get 

Eg{ug) = Eg'-{ug) + Eg^{ug) P Eg^iug). 

Thus it is enough to prove that for any k 

lim Eg^{ug) > ^ [ |(i/„(a;),? 7 fc)| d'H^{x). (3.10) 

V3 
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As the lattices and so the energies appearing in both sides of (3.101 are 
invariant under rotations of tt/S we confine our attention to prove (3.101 
for k = 1 and rji = ei. 


Note that the functionals consist of a superposition of 1-dimens- 
ional discrete energies related to the sublattices eZ x {m^/S/2} for m 
even varying in Z and to the lattices £(Z + (1/2)) x {my/S/2} for m odd. 
Hence, as a first step we will use the usual separation of scales argument 
on these 1-dimensional discrete energies in order to rewrite them in a 
suitable 1-dimensional integral form and then we glue the information 
back to the 2-dimensional setting by a slicing procedure. More precisely, 
for any m G Z set 


5^ = K X [m^/2, (m -b 1)^/2). 


Note that the stripes Sm n fl give a partition of ff. For any s G (0,1) we 
will construct a sequence {w|} C 5i?F(n,R^) with 1 dimensional profile 
along the direction ei and depending only in the first variable in each 
stripe. Thus let s G (0,1) be fixed and set 


/s = {i G Te n n : J(|u(i) — u(i -I- £ei)|/£) > sJoo} 


(for the sake of notation we drop the dependence on s in what follows). 
For any m G Z define We on R x {m^/S/2} to be equal to the value 
on (i, i -|-£) X {m^/3/2} if i G /e and to be the affine interpolation 
of the values Ui;{±),u^{± + £ei) on any other interval (i, i-|-£) x {m\/3/2}. 
Eventually, extend We on Sm n as We{xi,X 2 ) = Ws(xi, m\/3/2). 
Arguing as in the proof of Proposition |3.1| we infer that 

/ |V«;e|"dx-b£#(/E) <c (3.11) 

Jq 

and that Wj —>■ u in L^(fl,R^) and, up to subsequences, also almost ev¬ 
erywhere in fl. In addition 

Ee^{Us) > sJoa£#{Ie)- (3.12) 


Taking into account that 

#( 7 e )= Y.#{E(^Sm) 

m^Z 


and also that 

#{hnSm) = #{S{wt^'^)) 
for all y G Sm n , (|3.12 1 can be rewritten as 


2 

Es (r^e) ^ sJoo ^ 

v3 


f #{S{wt^’^))dH\y). 

in'i 


(3.13) 


By passing to the liminf in both sides of ( |3.13[ ) and by applying Fatou’s 
Lemma we get 


lim ^^^(me) > sJoo-^ [ liminf #(5'(w/^’^)) d77^(i/). (3.14) 

£^■0+ V3 Jn=i £^o+ 


On the other hand, estimate ( |3.11 1 ensures that for "H^-almost every y G 
the sequence We'^’^ is precompact in 5'i3’F(n®'^’^,R^) and converges 
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in measure to Hence by the 1-dimensional analogue of the lower 

semicontinuity Theorem [m we infer that 


lim Es^{ue) > sJoo-^ [ #{S{u'^’'^))d'H^{y). 

e-»o+ V3 Jn=i 

Moreover, thanks to Theorem |1.6[ a) with <; = 1, we have 

f [ \{v^{x),e,)\dH\x). 

in=i J s„ 


(3.15) 


(3.16) 


Letting s —>■ 1 concludes the proof of claim ( |3.9| 
yields (3.8l. 


Eventually, Lemma |3.4| 
□ 


4 Upper estimates for a class of “opening 
cracks” 


In this section we show the optimality of the bound (3.8 1 provided in 


Proposition |3.5| for a class of functions with ‘opening cracks’ with respect 
to the anisotropies inherited by the lattice. 

In the next section more complex geometries will be taken into account 
and the occurrence of different phenomena will be highlighted. Actually, 
the representation of the limit energy seems to take into account several 
factors, not all of local nature. 

In order to clarify the parameters playing a role in the asymptotic 
behaviour of , we prefer to deal with the case where only two rotations 

, R~ are involved in the target deformation u first. In this setting the 
crack Su splits H into two regions , in general not connected, each 

one underlying a rigid motion. We will show that in this case the P- 
limsup of Ee is finite even if the request of orientation-preserving recovery 
sequences affects substantially the form of the limit energy. 

In this process a relevant condition that translates the positive-deter¬ 
minant constraint through the crack is the following: 


{u'^{x) — u {x),R^i/)>0. 


(4.1) 


In fact, the inequality above ensures that any small triangle crossing the 
fracture site maintains positive (non negative) area in the codomain. Ac¬ 


tually, due to the discrete environment, condition (4.11 has to be assumed 


for u G D, since the only triangles entering in the construction have sides 
parallel to the directions lying in S. 

We start by proving the following characterization of the P-limsup 
whenever Su is contained in a straight line and ( |4.1[ ) holds. 

Proposition 4.1. For G 50(2) and x,u,q^ G with u G D, let u 
be defined as 


u{x) = 


R^ X + if {x — x,v) > 0, X G Fl 

R~ X + q~ if {x — X, o} < 0, x G FI. 


(4.2) 


Assume that R'^ 7 ^ —R~ and for R} -almost every x G Su R^, q^, v satisfy 
the condition 

{u'^ {x) — u~ {x), R^u) > Q (4.3) 

then 

r-limsupE,(M) < / y}{v)dR} =2 J^R}{Su). (4.4) 

£-> 0 + J Su 
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Proof. We claim that thanks to the hypothesis R'^ 7 ^ —R we may assume 
the stronger separation hypothesis 

{u'^ (x) — u~ (x), R^v) > 5 > 0 (4-5) 


for some fixed 5 > 0 and for ^^-almost every x £ S„. 

Indeed, if (4.51 does not hold, it is enough to choose a vector v such 
that {v,R^v) > 0 (set for instance v = R'^u + R~v) and consider the 
sequence us defined replacing g"*" with + du in (4.2 L Clearly us ^ u in 
and the left-hand side term of (4.41 is continuous along such a 
sequence. Condition (4.51 plays a key role in showing that the sequence 
(uj) defined as the pointwise interpolation of u in the nodes of the lattice 
nfl is admissible and accounts for the desired F-limsup estimate. 
Thus, set Ue(i) = M(i) for any i € L^nfl. Clearly —>■ m in R^); 

it remains to check that the positive-determinant constraint is satisfied in 
any triangle contained in fl. 



Figure 3: (a) Triangles deformed by the pointwise interpolation 


Before proceeding with the computation we label the triangles ‘cross¬ 
ing the fracture site’. By rotational and translational invariance we may 
assume v = (0,1) so that = { 2 ; £ : x = * 2 }. Let W := [ 2 a: 2 /\/ 3 e] 

and set 

Ve = {m G Z : {me, %/3W/2) £ Le n 11 } 

and 

im = {me, y/iNse/2) for m £ T>s. 

Note that u coincides with a ‘positive’ rotation on the vertices of those 
triangles not contained in the strip 

5 = {x : V3N^e/2 < X 2 < V3(W + l)e/2}. 


Hence, in order to ensure that £ Ade{Q), it suffices to check condition 
(4.31 for triangles T with vertices respectively im, im im +eri^ and 

im, im+e??^, im + £(r?^ —(see Fig. 2). This leads to prove the following 
inequalities: 


{{u {im + erj^)-u {im))^,u~''{im + erj^) - u (im)) > 0, 
((u+(im + ev^) - M’'"(im + £(J7^ - ?7^))^,M+(im - u~{im)) > 0 

for m G'Ds. These can be compactly rewritten as 

u+(im + £??^) - u~{im)) > 0. (4.6) 
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We claim that these conditions are fulfilled for e small enough. Indeed, 
for any infinite collection of indices {mg} C Z with € T>e, up to 
subsequences, we may assume that —>■ ® as e —>■ 0+. As a consequence 




-(x) and —>■ with £2 = X 2 - Thus, 


U^{T. 

assuming that (4.6 1 are violated for such a sequence of indices will 
lead to a contradiction to (4.5 I when passing to the limit as e —>■ 0+. 

Eventually we compute the asymptotic value of Es{us) to show that 
(uj) is a recovery sequence for the surface energy (p{uu) dV}. Clearly, 
the energy contribution in Es{u^) reduces to those pairs of indices one of 
which of type im with m £ (Dg; i.e.. 


+ (i„)| 


) 


E,{u,) = ^ ej(l 

m^'D^ 

+ ^ + -y^)) 


m^'Dg 

2 'y ^ s{joo + 0 ( 1 )) 

2£(/oo^{?xi ^ : ijTi, 


^ + 0(1). 


(4.7) 


Passing to the limit as e —>■ 0 and plugging the equality 


lim #{m £ De : i™ £ Su} 

E-J-0 + 


in (4.71 the conclusion is achieved. 


n\s„) 


□ 


The next proposition generalizes the result of Proposition |4.1| in the 
case when only two rotations are involved in the deformation and Su 
consists of segments having normals in D. 

Proposition 4.2. Let £ 50(2) and be given with yf 

—R~. Let S be a connected union of segments with normal belonging to 
D and let 


u{x) = {R^x + q^)Xsi+{x) + (R x + q )xsi-{x) 

where is a Caccioppoli partition subordinated to S. Assume that 

for -almost every x £ S it holds 

{u^(x) — u~ (x), R'^i'{x)) > 0, (4.8) 

then 

r-limsup£;,,(u) < [ g}{y[x))dH^ =2 JocH'[Su). (4.9) 

e->0+ Js 


Proof. Up to an approximation argument we may assume that S = UfcgiS'fe 
with Sk = for £ Le n fl, having constant normal Uk lying in 

D and satisfying = i’^^^ for any k. Thanks to the connectness hypoth¬ 
esis and to condition (4.8 1 , arguing by perturbation with infinitesimal 
translations, condition (4.8 1 can be assumed to hold in the stronger form 
{u^{x) — u~(x), R^y{x)) > S, for a positive 5. We claim that a recov¬ 
ery sequence is provided defining UE(i) = ^^^(i) for i £ n 14. Indeed, 
for any triangle T £ Tf intersecting S there exists a unique k £ X such 
that either T n 5 is a single point in Ss, or T n 5 = [i, i ± eUk] for some 
i € I/£ n 14. In both cases one may perform the same computation as 
in Proposition |4.1| and deduce that Ue satisfies the positive-determinant 
constraint. Eventually a direct computation shows that lim,j^g+ Eg{u^) = 
J°oY^kex‘Pi’^k)R^{Sk). □ 
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The previous approach can be pushed further to obtain the optimality 
of the bound (3.8 1 also for Su consisting of a line with normal v In 
this case we need to impose the opening crack condition (14.31) in a stronger 


sense; i.e., (4.3) must be satisfied also along the two directions i'i,U 2 G D 
generating i/ in one of the simplex of the Wulff shape {ifi < 1}. 

Proposition 4.3. For G 50(2) and x, v, G let u be defined as 


u(x) = 


R'^x + if {x — x,v) > 0, X £ Q. 

R~ X + q~ if {x — X, o} < 0, x G fl. 


(4.10) 


Let ui,V 2 G D fee such that vlip{v) = Xoi + (1 — X)i '2 for X G (0,1)- 
Assume that for R}-almost every x G Su R^,q^,v satisfy the conditions 


{u'^{x) — u {x),R^Vi)>S 
for i = 1, 2 and for some fe > 0. Then 

r-limsupi?E(u) < / ip{v) dR}. 
£->0+ J S„ 


(4.11) 


(4.12) 


Proof. The claim will be proved by an approximation argument, exploit¬ 
ing Proposition |4.2| and the lower semicontinuity of the T-limsup. It is not 
restrictive to assume Su connected and with endpoints on dLl. Note that 
hypothesis (4.11) ensures that R'^ 7 ^ —R~. Let h > 0 be a hxed small 
parameter and select {5^}^ a connected union of segments with normal 
equal to or to V 2 such that 

lim [ ipfuhix)) dR} = [ (p{y)dR} (4-13) 

+ JSt, jS„ 


and limji^+oo dniSh, Su) = 0. By simply extending the polygonals above 
we may also assume that each Sh has endpoints on dit and splits it in 
two components such that 

lim \A'^A{x G n : {x — x^v) > Q}\ — 0. 

h —>- + oo 

For h G N let Uh be the piecewise rigid function defined as 
Uhix) = (R'^x q^)XA+ + {R~x q~)XA- ■ 

h h 

It is easily checked that uj, —>■ u in L^(fl;R^) and, by continuity, Uh 
satisfies hypothesis (4.11) with a smaller positive 5 and Su replaced by 
= Sh- Applying Propositionto each Uh we get 


P- 


limsupi5e(uh) < / ip{yh)dR^. 

£—^ 0 ”^ J Suu 


The conclusion follows taking property (4.131 into account. 


□ 


In the following we consider the case in which the jump set Su consists 
in a triple point with coordinate normals. 

Proposition 4.4. Let xo G LI be fixed and set 

fli := {a: G n : (x - xo, ?7i^) > 0, (x - xo, 772 ^) > 0}, 

LI 2 ■= {x G LI : (x — Xo, t] 2 ^) < 0, (x — xo, 773 ^) < 0}, 
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^3 := {* € n : {x — xo,? 7 i^) < 0 , (x — xo,773^) > 0 }. 

Let u be defined as u(x) = (*) u^{x) = -R*^x + for 

suitable R^,R^,R^ G 50 ( 2 ) and q^,q^,q^ G R^. Assume that u satisfies 
the ‘opening crack’ conditions along Su = Ufe(^^ 7 s O i 2 ); 

{u^{x) — u^{x), R^rji^} > 0, {u^{x) — u^{x),R^pi^) > 0, on Oi n 523 

(u^(x) — u^(x), i?^» 72 ^) > 0, {u^{x) — u^{x),R^P 2 ^) > 0, on 52 1 n 522 

{u^{x) — u^{x),R^ri3^} > 0, {u^{x) — u^{x),R^ri3^} > 0, on 523 n 522 

(4.14) 

and that u satisfies the further compatibility condition 

(u^(xo) - M®(xo), (n^(xo) - M®(xo))^) > 0 (4.15) 


on the triple point xq. Then 

r-limsupi?£(M) < / ip{v) dR}. 

£-> 0 + 


(4.16) 


Proof. Up to composing u with infinitesimal translations we may assume 
that, for any e > 0 small enough, the points xo —e/ 3 (j 7 i+ 772 ), xo+e/ 3 (r 73 + 
7?i),xo + 5 / 3(772 — 773 ) lie on the lattice n 52. Set 7 i£(i) = 7 r(i) for any 
i € I/E n 52. Clearly —>■ u in L^(52,R^); it remains to check that the 
positive-determinant constraint is satisfied in any triangle intersecting Su- 


Thanks to hypothesis (4.151 we have that there exists 5 > 0 such that all 


the scalar products in (4.141 are greater than 6. By performing the same 


computation as in the proof of Proposition |4.1| we may deduce the positive- 
determinant constraint to be satisfied on any triangle intersecting Su and 
not containing xq. The only triangle left aside is then the one having 
vertices xo -|- £/2\/3(772 - 771), xo -I- e/2\/3(?73 - 771), xo -I- £/2v^(772 -I- 773). 
A direct computation shows that 

(7i^(X0 -I- 5/3(772 - 773)) - U^(xo - 5/3(771 -I- 772)), 

(w^(xo -I- 5 / 3(771 -I- 773 )) - ^^(xo - 5 / 3(771 -I- 772 )))^) = 

(u^(xo) - u^(xo), (u^(xo) - u^(xo))^) -1- 0 ( 1 ). 


Hence the conclusion follows by hypothesis (4.151. 


□ 


As a conclusion of this section we note that we have proved that the 
T-limit of Re is described by the anisotropic fracture energy 


R(u) = f ip{n)d'H^ 
•tSu 


(4.17) 


on all u which are piecewise rigid deformations such that Su consists of 
a hnite number of lines meeting at triple points and the opening-crack 
conditions in Proposition |4.4| are satished. This description extends by 
continuity to all piecewise rigid deformations u that can be approximated 
in energy by sequences of piecewise rigid deformations Uh satisfying such 
conditions. 
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Figure 4: A layer corresponding to fracture in the reference configuration, and 
the corresponding admissibility angle for u"*" — u~ 


5 Macroscopical failure of the positive- 
determinant constraint 

In the previous sections we have computed the F-limit of our discrete sys¬ 
tem when the macroscopic configuration satisfies compatibility conditions 
on the fracture site and at meeting points of three (or more) fracture 
sites. Those conditions can be regarded as a positive-determinant con¬ 
straint on the fracture. In this section we will see how all those conditions 
can be removed. In this way macroscopic configurations “with negative- 
determinant fracture” can be obtained from atomistic configurations sat¬ 
isfying a microscopic positive-determinant constraint, at the expense of a 
strictly greater energy. 



Figure 5: Surface relaxation allowing for {u'^ — u (deformed configura¬ 

tion) 


5.1 Removal of compatibility conditions on the 
fracture site — surface relaxation 

The surface energy ( |4.17| l has been obtained under a condition that implies 
that the fracture can be achieved by highly deforming only a single layer of 
positive-determinant triangles, keeping all other triangles at the minimum 
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Figure 6: Multiple microscopic fracture 


of the energy. We now show how, still using a single layer of triangles to 


create the fracture, we may relax conditions (4.111 in the case when v ^ D. 


We only illustrate the surface relaxation with an example. Consider a 
linear crack surface with normal v — (1,0), which does not belong to D. 
In Fig. 1^ we highli ght a layer of triangles that may be used to obtain a 
fracture satisfying (4.111. If R'^ = R~ = Id then the possible values for 
g"*" — q~ lie in the convex angle pictured on the right of the figure. We may 
obtain any g'*" — q~ with {q~^ — q~ ,v) > 0 if we also deform the triangles 
neighbouring the ‘fracture layer’. A possible deformation is given in Fig.[^ 
where each such neighbouring triangle has one side compressed by a factor 
^/3/2 and one side compressed by a factor 1/2. This construction gives 
an estimate for the energy density on the fracture with 


2 + 2j(^)+j(‘) 


(while ipiy) = 2). 

Note that such a wide relaxation is not necessary for fixed q^, q~ if 
(g^ — q~ ,v) > 0. In general then, even in this simple geometry, we will 
obtain a surface energy density depending on — u~. 



Figure 7: Multiple microscopic fracture - reference and deformed configurations 
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5.2 Removal of compatibility conditions on the 
fracture site — multiple microfracture 

By introducing more layers we may remove all constraints on Su, while 
the macroscopic energy varies by a factor proportional to the number of 
additional layers. 

In Figure the macroscopic deformation on the fracture does not 
satisfy the conditions 

— u~ (x), Vi) > 0 

for the two “microscopic” coordinate normals Vi. However, the discrete 
functions in the construction depicted in Figurej^all have positive-determinant 
interpolations thanks to the introduction of a “fictitious” layer of atoms. 

In this case the limit energy per length doubles. 

A variation of this example in given in Figures and to highlight 
that the introduction of an extra layer of atomic interactions allows to 
remove the condition {u^{x) — u~{x),R^v) > 0 even when is a coordi¬ 
nate normal. By repeating this process the macroscopic deformation may 
exhibit a interpenetration phenomenon. 



Figure 8: Multiple microscopic fracture - macroscopic configurations 


Note that for some deformations surface relaxation may be energet¬ 
ically convenient with respect to mult iple-layer fracture. This clearly is 
the case of the example in Section 5.1 when — q~ is close to satisfying 


the compatibility conditions (|4.11[|. 


5.3 Removal on compatibility conditions on “triple 
points” — micro-deformed fracture 

The condition on triple points in the previous section ensures that the 
deformation of a microscopic triangle at that point is of positive de¬ 
terminant (and hence, being a single triangle, gives a negligible energy 
contribution). If such a condition does not hold then the use of (small 
variations of) pointwise interpolations on the different regions of the un¬ 
derlying partition is not possible, since for e small there will always be 
a microscopic triangle whose vertices are mapped in three points which 
fail the positive-determinant constraint. However, it is possible to use a 
different interpolation by introducing an additional microscopic fracture 
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Figure 9: Triple point with “negative determinant” - reference and deformed 
macroscopic deformation 


enclosing small sets where the deformation is not in SO{2). Note that 
this is possible if such sets have the dimension of an interface. 

In Figurej^it is represented a deformation with a triple point failing the 
positive-determinant condition. In this case we introduce a microscopic 
approximation as represented in Figure 



Figure 10: Triple point with “negative determinant” - reference and deformed 
microscopic deformations 


On the fracture site we introduce a segment [A, X] where the point- 
wise single-layer interpolation of the jump is substituted by a double-layer 
approximation. Note that the image [A', X'] of this segment undergoes an 
additional linear deformation. Note that the energy of this approximation 
provides the additional contribution 

where the second term is due to the compression of the triangles on the 
segment [A, X] in the reference configuration. This energy depends on 
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the choice of X and X', which are variables in the construction (sat¬ 
isfying some constraints due to the positive-determinant requirement in 
the resulting construction). This shows that even in this simple case an 
optimization problems arises between the introduction of an additional 
microfracture and a microscopic compression. Of course, more complex 
constructions with more parameters can be also introduced. 

5.4 Global failure of impenetrability constraints — 
optimal decomposition and healing microfractures 

In the constructions illustrated above we could exhibit a microscopic re¬ 
covery sequence by working separately on each fracture site or triple 
points. In the presence of a complex geometry of the domain, besides 
a use of those constructions one has also to take into account the possibil¬ 
ity of introducing further “fictitious” microscopic interfaces to get around 
impenetrability constraints. As a simple example, we may consider the 
deformation in Figure m where the central smaller triangle is removed 



Figure 11: Deformation violating the impenetrability constraint 


and translated from its position in the larger triangle. This macroscopic 
deformation can be approximated by microscopic ones all satisfying the 
positive-determinant constraint (see Figure [I^ . One such approximation 
can be obtained by translating a rhombus and subsequently composing 
this translation with a deformation rotating half of this rhombus as de¬ 
scribed in Section 5.2 (see also Figure]^. Note that this last rotation 
entails the introduction of one or more microscopic fictitious layers of 
atoms. 

Another simple example is depicted in Figure [T^ where the triangle 
to be removed is strictly contained in the interior of a larger triangle. In 
this case, in order to proceed as in the previous construction, one has 
to introduce an “auxiliary fracture” as the segment [X, X'] in the figure. 
The determination of the optimal shape and location of such auxiliary 
fracture sites is clearly a complex optimization problem. Note that in this 
case, also the determinant constraint for triple points has to be taken into 
account. 
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Figure 12: Construction of approximations satisfying the positive-determinant 
constraint 



X’ 


Figure 13: Deformation that can be obtained with auxiliary fractures 
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6 Necessary conditions for opening cracks 


In Section we have shown a lower bound with an anisotropic Griffith 
fracture energy, which is optimal on a family of displacements with an 
opening-crack condition on the fracture site (Section . We now show 
conversely that if the limit energy at a point xo € Su is not greater than 
the lower bound then necessarily the function u satisfies a opening-crack 
condition. 

We now consider discontinuity points where the fracture energy density 
is minimal; i.e., the inequality in (3.8 1 is sharp, and derive necessary 
conditions on the crack opening. To this end we introduce the measures 


^ ( 6 . 1 ) 

(i.DeATan) ^ 


These measures are nothing else than a way to measure locally the energy 
Fs{Ue). If fJ-e is a bounded sequence of measures and sup^ ||'ii£||Loo{n,E 2 ) < 
-too, then, arguing as in the proof of Proposition |3.1[ we infer that is 
precompact in and that —>■ u in We also suppose that 

the weak* limit jj. of /ig exists. 

Proposition 6.1 (necessity of an opening-crack condition). Let Ue —>■ u, 
and let xo G Su be such that 


d/a 

dW L Su 


(xo) < (p(nu(xo)). 


( 6 . 2 ) 


Then we have 

(u'^ixo) - u~(xo), R'^i'uixo)} > 0, (6.3) 

where G SO{2) are the two constant matrices coinciding with Vu on 
both sides of Su at xq. 


Proof. For the sake of brevity we will denote vq = i’u{xq). 

Let be a square centered in 0, with side length 1 and an edge 
orthogonal to vq. With fixed p > 0 and y G <3''°, let 


= Ue{xo+ py), 


which, by definition of Su converges as e —> 0 and p —> 0 to the function 

-/ N I u+{xo) if {y,vo) > 0 

u(y) = < 

' u (a:o) if {y,vo) < 0. 


By the weak* convergence of and (6.21 we deduce that 


Ee{Ue,Q''p°{xo)) < PT{^) + Op(l) + 0^(1) 


where Q''p°{xo) is a cube centered in xo, with side length p and an edge 
orthogonal to vq, and Op(l) and 0^(1) are infinitesimal as p —>■ 0 and e —>■ 0, 
respectively. 

We fix s > 0 and set 

Sp = G : I I > g Iqj. least two sides (i, j) of r|. 

We claim that we may connect the two opposite sides of Qp° (xo) parallel 
to Vo with a path {R : i = 1,... ,M} (depending on p and e, but we 
omit such a dependence for the sake of notational simplicity) consisting 


27 






of triangles such that Ti and T^+i have a common side and Ti € up to 
a number of indices that is o{ple). Indeed, note that for any T ^ SI we 
have 

I I ^ 2s for every side (i, j) of T. (6.4) 

If no path as above exists then we may construct cpje disjoint paths in 
Le n Q''p°{xq)-, i.e., sets of indices {i^ : n = 0,..., Mj} with and 

€ eS, such that 

(ij -xo,vo) < -^p + £, {±\j. -xo,i'o) >^p-£ 


and 




< 2s for every n 


(6.5) 


(for a construction of such paths we refer to the proof of Theorem 4(ii) in 
|10p. Since — u^{xo) converge to 0 close two opposite sides of Q''° it 
is not restrictive to suppose that 

^e|u'"(i^) -u“(2;o)| + ^£|w^(iifP - u+(a;o)| =po^(l). 


By (6.51 we have 


cp|u+(a:o)-w (a:o)| < ^ e(|w’^(a;o) - u'"(iMp| 


+|w^(iMj) - w^(io)l + l«^(io) - “ (®o)|) 

Mj 

< pOe(l) + ^^e|u‘"(ii)-u^(i^_i)| 

3 ^=1 

Mj 

< pO.(l)+^^£^ 2s < poe(l) + p^2s. (6.6) 

j n=l 

Dividing by p and letting e and p tend to 0 we then obtain |M"'"(a;o) — 
u“(a;o)| = 0, which gives a contradiction. 

We consider the connected set 

M 


Note that by the convergence Vs,p — >■ u, for all 5 > 0 the set Tg n Q'3‘{xq) 
will be contained in the strip {x ■. \{x — xq,vo)\ < 5p} for e sufficiently 
small. It is not restrictive to suppose that Q'3’{xo) n is composed 
exactly of two polygonal chains Pt and P+. Note that the set 

{Ti : dTi n PI 7 ^ 0} 

still gives a path with the same properties of {Ti}. We can therefore 
suppose that this is our original path. 

The contribution of the energy restricted to interactions in gives 


E,{u,,Q';%xo)r)Te) > J{s)<p{uo)p{l-5) 


(6.7) 


(see Proposition 3.51. We then deduce that 


#{r £ SI-.Tit TJ < c^5 


( 6 . 8 ) 
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and, in addition, there exist disjoint polygonal chains joining PL and 
the side of Qp°{xo) lying on {y ’■ {y — xo,v) = —pji} such that they do 
not intersect S^. Since Vs,p —>■ u, arguing as in ( |6.6| l, we deduce that 

[ \Ue - u'^ixo)\d'H^ = pOs(l) + o{p) (6.9) 

'>p± 

as £ —>■ 0. 

For Hxed p > 0 we may now consider the scaled functions 

wi{y) = -Ue{xo + py) = -U£.p(y). 

P P 

We may still argue as in Proposition |3.1| and deduce that converge to 
u'‘{y) = ^u(xo + py). Note that Vu’’{y) = Vu{xo + py). 

For Hxed p, we may now consider the functions 


We = w^Xq-ovIsi’ 

which are now discontinuous at dSl- They satisfy 

(i) the determinant of VwC is equibounded, by the definition of 

(ii) VWj are equibounded; 

(hi) l-L^idSL) is bounded; 

(iv) the measures L dS" are equibounded, by (ii) and since 

I I < 2s by definition. Here and below we simply denote by u the 
normal to a discontinuity set without specifying of which function, which 
is clear from the context. 

By Theorem 5.8 in [3] from (i)-(iv) above and the convergence We —>■ 
we deduce that the measures L dS^ converge to L SuP ■ 

We choose the orientation of so that the determinant constraint on Ue 
can be rewritten on each segment of P± as 


where uj and are the values on two vertices of the corresponding 
triangle, with G Pt and ut G P+- We then have, by (6.91 

dUe 


0 < 


jQ7(-o)nPt\9uP / 


^ L 


(u+(a:o) - u {xo))-^\d'H^ + poe(l) + o(p) 

p [ (^^,{u'^{xo)-u~{xo))^)d'H^ +pOe{l) + o{p) 

jQ'^onj(Pi-xo)' 

P [ (^■'"(a^o) - w“(a:o))^)dH^ + pos(l) + o(p) 


/Q‘'onS„p 

P {R~ , {u'^ (xo) - U~{xo))^) + pOe(l) + o(p) 

p {R~IP, u'*'{xq) - U~{xo)) + pOe(l) + o(p). 


Dividing by p and letting £, p —> 0 we obtain the first claim in (6.3l. The 
analogous inequality with R^ip is obtained in the same way. □ 


Remark 6.2 (conjecture). We have obtained a weaker necessary condi¬ 
tion than the sufficient ones in (4.111. We conjecture that this is only a 
technical issue and indeed those condition are also necessary in order that 
the energy density be minimal on Su- 
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7 Small deformations 


In the previous sections we have showed how the “positive-determinant” 
constraint on the fracture is not closed, and must be substituted by com¬ 
plex non-local energetic considerations. The constraint simplifies in the 
case of small deformations; i.e., if we require that 

u = id -I- SeV (7.1) 

(or any function with gradient in 50(2) in place of the identity) with 

6 « 1 . 

If « 5e then we can still apply the compactness argument as above 
to obtain, from the boundedness of Eelus), that v determines a partition 
into sets of finite perimeter. Since u'^{x) — u~{x) = 5^{v^ {x)— v~ (x)) and 
= id -I- o(l), in the case of a single coordinate normal u the constraint 
on u 

{u^{x) — u~ (x), R^v) > 0 
can be translated into one on v that reads 


(u^(a;) — V (x), ly) > 0 


(7.2) 


(in the case of a general non-coordinate normal we have to take into 
account the microscopic oscillations of Sv as usual). Note that this is 
the usual infinitesimal opening fracture constraint, which is a closed con¬ 
straint. 


Remark 7.1 (linearized regime). The case 5^ of the order of corre¬ 
sponds to the one studied by Braides, Lew and Ortiz in the one-dimensional 
case when bulk and surface terms have the same scaling. Hence, the 
compactness argument used above does not hold, and the limit energy 
possesses both a bulk and a surface term. The linearization of the bulk 
term has been studied by Braides, Solci and Vitali m, and can be used 
together with the analysis above. The resulting T-limit can be explicitly 
written as 


CVu : Vu dx -I- 


L 




on piecewise functions whose discontinuity set is a union of smooth 
lines, V satisfies the infinitesimal opening fracture constraint on 5„. Note 
that a construction as in Section [5.3| implies that the condition on triple 
points can be dropped altogether (since microfractures can be chosen van¬ 
ishing with 5g). The proof of the upper bound for the T-limit is a techni¬ 
cally complex matter due to the lack of density theorems for v satisfying 
the resulting constraints, and will not be addressed here. Related results 
can be found in works by Friedrich and Schmidt da mini HZ]- 


8 Conclusions 

We have examined a two-dimensional system of nearest-neighbour inter¬ 
actions for Lennard Jones pair interactions parameterized on a triangular 
lattice, with a microscopic positive-determinant constraint which mimics 
the effect of long-range interactions and limits the ground states to rota¬ 
tions. The goal is to exhibit a continuum approximation so as to test the 
validity of such a simplified model. 
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In parallel to the one-dimensional case we have focused our attention 
on the fracture term of the resulting continuum approximation by suitably 
scaling the energies. We have thus determined 

• a surface energy density ip which reflects the triangular symmetries 
of the underlying lattice; 

• conditions on the interface that reflect the positive-determinant con¬ 
straint. The conditions on the interfaces are of a novel type that take into 
account both the gradient of the deformation on both sides of the fracture 
and the orientation of the fracture site in the reference configuration. In 
addition we also have a positive-determinant constraint on points where 
more fracture meet (triple points); 

• such conditions are not a closed constraint, and can be removed 
by adding “hctitious” micro-fractures. The optimal location and form of 
those micro-fractures depends on the corresponding macroscopic defor¬ 
mation, and is a complex optimization problem; 

• in the case of small deformations the conditions on the interfaces 
reduce to the infinitesimal opening fracture constraint, which actually 
suggests a linearized approximation in the spirit of Braides, Lew and Ortiz 
with that constraint on the fracture (also analyzed in recent works by 
Friedrich and Schmidt [TS1II3I1S1II3)- 

We have neglected interactions other than nearest neighbours. In this 
way the interfacial energy does not reflect the possibility of surface re¬ 
laxation; i.e., the fact that atomistic interactions are unbalanced close 
to interfaces. This is an important effect, especially at the boundary of 
the domain and in the determination of the location of fractures. It has 
been partially addressed by Theil [^, and should be included in further 
investigations on the subject. The one-dimensional analysis by Braides 
and Solci m suggests that considering longer range of interactions (e.g., 
next-to-nearest neighbours) could be used in the place of the determinant 
constraint in order to eliminate non-local effects for non-opening cracks, 
at the expense of introducing internal and external boundary layers. 


Appendix: positive vs non-negative deter¬ 
minant constraint 


We have considered a strictly positive microscopic determinant constraint. 
Strictly positive inequalities are weakened in the limit; however, our choice 
of not directly considering weak inequalities allows to rule out some ad¬ 
ditional “unphysical” deformations which would have to be taken into 
account by directly considering a non-negative microscopic determinant 
constraint. This would correspond to allowing microscopic interpolations 
to “collapse” triangles to segments on the jump set even though such 
collapsed triangles cannot be viewed as a limit for the strictly positive 
determinant case. This may happen in the case of rotations £ 5*0(2) 


on both sides of the interface with = —R . In Fig. 14 we depict such a 


microscopic deformation along a single coordinate line, where all points of 
two rows are aligned. Note that allowing zero-determinant deformations 
would include such a macroscopic deformation in the set of “minimal in¬ 
terfacial energy”, while in our setting the same must be achieved by the 
introduction of at least one extra layer of atoms, thus doubling the energy. 

Another possibility for having = —R~ is with the partition com¬ 
posed of a pair of supplementary angles as in Figure |15| In this case it 
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Figure 14: Zero-determinant fracture with a planar interface and a 180-degree 
rotation -- reference and deformed microscopic configurations 


is possible to “rotate” one of the two angles by tt with positive or zero 
determinant on the deformed triangles along the discontinuity set. In 
both cases, however, we have a macroscopical failure of impenetrability, 
so we have regarded these cases as degenerate by considering only the 
strictly-positive-determinant constraint. 




Figure 15: Zero-determinant fracture with an angular interface and a 180-degree 
rotation - reference and deformed microscopic configurations 
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